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This paper briefly reviews the concatenated Reed-Solomon (RS)/Viterbi coding 
system. Then the performance of the system is analyzed and results are derived with a 
new simple approach. We present a functional model for the input RS symbol error 
probability ir. Based on this new functional model, we compute the performance of a 
concatenated system in terms of RS word error probability, output RS symbol error 
probability, bit error probability due to decoding failure, and bit error probability due to 
decoding error. Finally we analyze the effects of the noisy carrier reference and the slow 
fading on the system performance. 


I. Introduction 

The Voyager 2 spacecraft, which was launched in 1977 and 
has encountered Jupiter and Saturn systems in 1979 and 1981 
respectively, will fly by Uranus in 1986 and Neptune in 1989. 
Voyager delivered 115.2 kbps at Jupiter (5 AU from earth), 
and 44.4 kbps at Saturn (10 AU from earth), both with 5 X 
10“^ bit error probability. In order to enhance Voyager’s 
communications capability, NASA’s Jet Propulsion Labora- 
tory is expected to exercise a plarmed option by switching on 
its Reed-Solomon (RS) encoder on the Voyager spacecraft. 
This results in a RS/Viterbi concatenated coded communica- 
tion link to increase the achievable data rate. Reed-Solomon 
decoders will be installed in NASA’s Deep Space Network 
(DSN) ground stations in time for the 1986 Uranus encounter. 


The concatenated RS/Viterbi channel coding and its perfor- 
mance has been considered in the past; for example, see Refs. 
1-5. This paper gives a brief review of this topic and shows the 
derivation of the results with a new simple approach that will 
be useful to system design engineers. 

In concatenated RS/Viterbi channel coding, the key param- 
eter for evaluating the performances of the system is the input 
RS symbol error probabihty n. Based on previous simulation 
and measurement results, we give a simple functional model 
for RS decoder input symbol error probability n. Using this 
model, we evaluate the RS word error probability, output RS 
symbol error probability and bit error probability. Finally we 
consider the effects of phase jitter and slow fading channels: 
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namely, the effects of Rayleigh, Rician, and log-normal fading 
channels on the performance of this concatenated RS/Viterbi 
channel. We use the outer code RS (255, 223), convolutional 
rate 1/2 and constraint length 7 code to illustrate these results. 


II. System Model 

The block diagram of the concatenated coding system is 
given in Fig. 1. The Reed-Solomon/Viterbi concatenated code 
consists of a Reed-Solomon (RS) outer code and a convolu- 
tional iimer code which is Viterbi decoded. 

The binary input data sequence is divided into J bit 
sequences to form symbols over a fl/=2'^-ary alphabet; i.e., 
there are M possible RS symbols. The RS coder then encodes 
the symbols such that any combination of t or fewer symbol 
errors per RS word (2-^- 1 symbols per word) can be cor- 
rected. A very simple block diagram of an RS block coder is 
shown in Fig. 2. All coding and decoding operations involve 
RS symbols, not individual bits. Here K = M - \ - 2t infor- 
mation symbols (or J[M - 1 - 2r] information bits) from 
some data source enter the RS coder to the left. The result of 
coding operations is a codeword of length A^= 2-^- \=M- 1 
symbols, of which the first K = M - 1 - 2f are the same 
symbols as those entering to the left. This makes the code 
systematic. The remainder of the codeword is filled in with 2t 
parity symbols; t= (N - K)/2 represents the number of cor- 
rectable RS symbol errors iri an RS codeword. That is, if t or 
less RS symbols are in error in any way, the decoder will be 
capable of correcting them. An RS symbol is in error if any of 
the J bits making up the symbol are in error. The minimum 
distance of RS code is 2r + 1 symbols. ' 

The interleaving buffers are required because the inner 
Viterbi decoder errors tend to occur in bursts, which occasion- 
ally are as long as several constraint lengths. Without interleav- 
ing, Viterbi decoder burst error events would tend to occur 
within one RS codeword. Thus over a period of time there 
would be a tendency for some codewords to have “too many” 
errors to be corrected. The performance of the RS decoder is 
severely degraded by highly correlated errors among several 
successive symbols. The purpose of interleaving and 
de-interleaving is to make the RS symbol errors, at the input 
of the RS decoder, independent of each other and disperse the 
RS symbol errors, in other words, to break the burst errors out 
of the Viterbi decoder among several code words. 

The level of interleaving / corresponds to the number of RS 
code words involved in the interleaving and de-interleaving 
operation. Interleaving and de-interleaving operations over a 
Viterbi chaimel can be explained simply by considering two 
IX 2-^ - 1 matrices, one at the input of the channel and one at 


the output. For interleaving, put the code words with length 
2-^ - 1 in rows 1, 2, • • • , / of the matrix, then transmit the 
symbols of columns 1, 2, • • • , 2-^ - 1 through the chatmel. For 
de-interleaving, do the reverse operation. For the RS code 
(255, 223), simulation results (Ref. 3) have shown that the 
interleaving level of /= 16 is sufficient to make the RS symbol 
errors independent of each other for the (K = 7, r= 1/2) 
Viterbi channel. Thus we can assume henceforth that the 
combination of interleaving, convolutional code, additive 
white Gaussian noise (AWGN) channel, Viterbi decoder and 
de-interleaving creates an equivalent Af-ary discrete memory- 
less channel (DMC) with transition error probabihty tt, shown 
in Fig. 3. 

ill. A Functional Model for the RS Decoder 
Input Symbol Error Probability 

Based on simulation results of Linkabit for tt given in Ref. 3 
and then using a least-squares curve fit, we proposed a simple 
model for n as 



where /3 q = 4.9551, = 5.2275, T* = max (0o/^i> ^)> T is 

Viterbi decoder node synchronization threshold (for EJNq 
< T Viterbi decoder produces random output), and EJN^ 
is the bit SNR of the Viterbi channel (inner code). For perfect 
node synchronization T*= Figure 4 shows the baseline 

symbol error probability n for the (255, 223) RS code vs 
E^Nq of the Viterbi channel (inner code), using (1). 

IV. Error Events in RS Decoding 

Consider an {N, K) Reed-Solomon code. We know that 
r-error correcting RS code can decode a code word correctly if 
the number of symbol errors are t or less where t = (N- K)/2. 
Consider the space of all received vectors y of length N and 
spheres of radius t around code words in this A-dimensional 
vector space, as is shown in Fig. 5. These spheres do not 
overlap because the minimum distance between codewords is 
2f + 1 . Since we have a linear code, without loss of generahty, 
we assume the aU zero codeword Xq is sent. In decoding the 
RS code words, depending on what region in the observation 
space of Fig. 5 the received vector y falls, three disjoint events 
may occur. The first event occurs if there are t or less input RS 
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symbol errors in a received codeword. This occurs if the 
received vector y falls in the double-shaded sphere around x^, 
ye5^. In this case the decoder successfully corrects the errors, 
decodes Xq as the codeword and outputs the correct informa- 
tion block Uq = Uq. Then, the “correct decoding” event 
occurs. The second event happens if there are more than t 
input RS symbol errors in a received codeword and the cor- 
rupted codeword is not within a distance of t symbols to any 
other codeword. In other words, y does not fall in any sphere 
of radius t around the codewords, ye5^. In this case the RS 
decoder fails to decode and we say that a “decoding failure” 
event occurs. Here, the decoder just outputs the first N- 2t 
undecoded Viterbi channel symbols (input RS symbols) that 
may contain symbol errors. The third event happens if there 
are more than t input RS symbol errors in a received codeword 
and the corrupted codeword is within a distance of t symbols 
to some other codeword than the correct codeword x^. This 
happens if y falls in any one of the single-shaded spheres, yeSj. 
Here the decoder decodes incorrectly and outputs an incorrect 
information block. Then, a “decoding error” event occurs. 


1. Derivation of Pgp- Recall that when the RS decoder 
fails to decode, yeSp, the decoder simply outputs the first M - 
1 - 2r Viterbi channel symbols as the “decoded” information 
block Uq. Hence the output RS symbol error pattern is exactly 
the same as the first M- 2t undecoded input symbols to 
the RS decoder. An output RS symbol at position say k,k = 
1, 2, • • •, M- 1 - 2f is in error, when the input RS symbol at 
the same position k is in error and there are at least t input RS 
symbol errors in the remaining A^- 1 symbols of received 
codeword (i.e., other than the position k). The probability of 
the input RS symbol at position k being in error is tt. The 
probability that there will be t or more input RS symbol errors 
in the remaining positions, i.e., the A^- 1 places of the received 
codeword, is 

P^(;v - 1 , r - 1) = ■ ^ r ^ 1 TT-d - (4) 

m=t ' ^ 

Thus 


V. Error Performance Analysis 

In this section we derive expressions for RS codeword, RS 
symbol, RS bit, and RS information block error probabilities. 

A. RS Codeword Error Probability 

Denote the word error probability for an (N, K) RS code 
by t). An RS codeword is in error when there are more 

than t chaimel symbol errors in a received codeword y; i.e., 
ye {Sp U Sj} . Therefore, we should consider all combinations 
of k channel symbol errors within the N symbols, for all k > t. 
Thus 


Ps,F = ( 5 ) 

2. Derivation of P^j. As we have seen, a decoding error 
occurs if the received codeword y belongs to one of the 
spheres in 5). Consider a particular codeword x„ which has 
exactly n particular nonzero symbols. This codeword has a 
distance n from the all zero codeword x^ and hence has 
weight n. Let 

W W W ■ • ■ W 

be n independent identicaUy distributed (iid) random variables 
corresponding to the n nonzero symbols of x„, and 


Pw(N,t) = 


so 

fc=f+i 


7r*(l - 


( 2 ) 


where tt is the RS decoder input symbol error probability 
leaving the Viterbi decoder (groups of J bits). 


B. RS Symbol Error Probability 

Denote the RS symbol error probability by Pg. This error 
may result when the received codeword is either yeSp or yeSj. 
Denote the RS symbol error probabihties when yeSp and yeSj 
and Pg j, respectively. Since these events are disjoint, 

we have 


P +P 


( 3 ) 


■ ■ ’^N-n 

he N - n iid random variables corresponding to the A^ - « zero 
symbols of x„. Each of these random variables W^■, i='l, 
2, • • • , n, and Z^;/ = 1, 2, • ■ • ,N- n, can be either 0 or 1. Let 
us compute the probability that Xq is sent and the decoder 
decodes it as x„, i.e., the decoder has made an incorrect 
decision. Denote this probability by Pt{Xq ^x„}. This event 
occurs if y falls in a sphere with radius t around the codeword 
x„. The probability that a zero symbol is sent and a nonzero 
symbol of a particular value is received is 7rl(M- 1). But the 
probability that a zero symbol is sent and a nonzero symbol of 
any value is received is tt, since there are Af- 1 nonzero 
symbols. For computation convenience, without loss of gen- 
erality, we can assume the first n symbols in x„ are the 
nonzero symbols. Now let y = (yj • • • ,yjy) and x„ = (jc, „ . 
^ 2 ,n . • • • . ^jv,n)> 2 nd define 
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0 ; = X, 


i,n 


W. = 
1 


1 ; j', 


i = 1, 2, • • • ,n 


hn 


then given Xq is sent 

= 01 = 

Similarly, if we define 


Z. = 

I 


® ’ ^n+i ^n^i,n 


1 ; y„, , ^ X 

' m-i n+i,n 


i = 1, 2, • • • ,N- n 


then given Xq is sent 


Pr{Z. = 0} = 1 - 7T 

Pr{Z^ = 1} = 7T 


- • - :n 

tt- y fT' . 

where 
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( 6 ) 


(7) 


( 8 ) 


(9) 


and 


Pr{Z=j}= a^(l (13) 


Therefore 

t mm(t-iN-n) 

E Pr{W = i}Pr{Z=j} {\A) 
1=0 /=0 

This is an exact analytic result. An upper bound using the 
Chernhoff bound technique (Ref. 6) is derived in Appendix A. 
This gives a simpler expression than the exact result in (14). 

Let q(n) denote the number of codewords of distance n 
from x^; then the symbol error probability Pgj due to decod- 
ing error is 


Let 


N 


(15) 


n=2f+l 


N-n 


w 


= 5^ and Z = ^ Z. 


/=t 


/=! 


RS code is a maximum distance separable (MDS) code. 
Therefore we have (Ref. 7) 


Then denoting the weight between two vectors by cj(*,*) we 
have 


q{n) = (M-\) (16) 

m=0 


= Pr{W + Z<t] 

But 

t 

Pr{W + Z<t} = /V{H/ = /}/V{Z<f-r} 

/=o 

t mm(t-i,N-n) 

= '^Pr{W=i] Y. Pr{Z = f} 
i~0 /=0 


C. RS Bit Error Probability 

Denote the RS bit error probability by fj,(/?i). This bit 
error is either due to decoding failure or decoding error. 
Denote the RS bit error probability when yeSp and yeSj by 
Pi^p{R!^ and P^^(RS), respectively. Since these events are 
disjoint, we have 

P,(RS) = P^piRS) + P^/RS) (17) 

( 11 ) 

1. Derivation of fj,p(/?S). Recall that when the RS 
decoder fails to decode, it outputs the first M- 1-2/ input 
RS symbols as “decoded” information block Uq. A bit at 
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position say m, m = \,2, • • ■ ,J within the RS output symbol 
at position say k\ k = 1,2, ■ ■ • ,M - l-2tofa codeword is in 
error if the bit at position m within the input RS symbol at 
position A: of an undecoded codeword is in error and there are 
/ or more RS symbol errors in the received codeword out of 
the Viterbi decoder in other positions than position k. But the 
probability that a bit would be in error in an undecoded RS 
codeword is the bit error rate of the Viterbi decoder denoted 
by (Viterbi). The probability that there are t or rhore RS 
symbol errors in an undecoded codeword in positions other 
than position k is given by (4). Thus the RS bit error rate 

Pi,p(RS) = (Viterbi) • P^(N- 1, t- 1) (18) 

2. Derivation of P^^ j(RS). 'Ne have derived the symbol 
error probabUity in (15). Now the bit error probability 

=2-#T)'’w 

Based on our functional model for ir, /(RS) is shown in 
Fig. 6 for (255, 223) RS code using (19) with (14) and (15). 
In Fig. 7 we have shown P^ p(RS), Pg p and P^ for concate- 
nated coding with (255, 223) RS code and (7, 1/2) convolu- 
tional code vs the concatenated channel. For 

Eb/^o 2dB, Pw(^~ 1. f ■ 1) - 1; hence (18) becomes 
Pj^p(RS) = Fj(Viterbi). Note that the code rate of the RS 
code is 

RS code rate = (N - 2t)IN (20) 

As is seen from Figs. 6 and 7 


Using this fact, we determine the RS of information block 
error probability. 

D. RS Information Block Error Probability 

Note that if a codeword is in error it is not necessary that 
the corresponding information block be in error. This error 
probability is important for source coding (data compression). 
Denote the information block error probability by P^N, t). 
Ignoring the probability of decoding error event to find 
Pj(M, t), we should subtract the probabUity of all possible 
patterns of symbol error that happen only in the parity check 
symbols, from Pw This imphes there is no error in the 
information block. Thus we have 

2t 

Pj(N, t) = P^(N, t) - ^ (1 - (23) 

Ac=t+1 

Obviously, in the practical range of error probability the last 
term can be ignored and we can say approximately 

Pj(N,t)^P^(N,t) (24) 

VI. Effects of Phase Jitter 

If the data rate is sufficiently large with respect to the PLL 
loop bandwidth so that the phase error does not vary signifi- 
cantly during the Viterbi decoder error bursts, then the phase 
error can be assumed to be constant during one RS symbol. In 
addition, since we have /= 16 level of interleaving, the phase 
errors affect the RS input symbols approximately indepen- 
dently. Noting the approach we have taken for the derivations 
of and P)^(RS), we have 


and 


P «P 

^S,I S,F 


P^j(RS)«P^p(RS) 


( 21 ) 


P^ = E{P^(<p^ . ^2. • • ■ ^;v)> = ^ 

*=f+i 

(25) 

and 


Therefore 


and 


P ^P 
^5 ~ S,F 


P^(RS) = P^(Viterbi\<l>)'^ r<0)'^ [1 - 

fc=f ' ' 

(26) 


P^(RS) - Pf,p(RS) (22) 

At this point, since the probability of a decoding error 
event is very small, it can be ignored for further computations. 


where 


rr(?y =1 ;r| 

•'-TT 


AC. 


cos^ 0 I p(<t>) d(j> 


(27) 
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and p(0) is the probability density function of phase jitter 
given by (Ref. 8) 


P(0) = 


exp (p cos (f) 
27r/„(p) ’ 


-7T<0<7r 


(28) 


where p is the SNR in the tracking loop bandwidth and is 
given b^), replacing EJN^ with (EjN^) cos^ <t>. Figure 8 
shows 7i(0) for various loop SNR p. The effect of a noisy 
carrier reference on Reed-Solomon Viterbi bit error rate per- 
formance, give n by (26), is shown in Fig. 9. Note that in (26) 
{Viterbi\<^ can be computed from 


where m{t) contains the telemetry information which can be a 
subcarrier, biphase modulated with a binary data stream. Now 
if this signal is passed through a fading charmel, at the output 
we have 


r{t) = y/2a(t)%m [co^r + /n(r) + 0(f)] + «(0 (32) 

where a(t) is random amplitude process and 6(f) is random 
phase process. If a{t) and 0(f) change slowly with time, and 
the spectral bandwidths of a{t) and 0(f) are narrow with 
respect to average loop bandwidth of PLL at the receiver, then 
0(f), together with the carrier phase, can be tracked by PLL. 


{Viterbi\<p) = J* cos^ p(<t>)d4> 


where 


Here we consider a slow varying fading channel with perfect 
tracking. Depending on the channel, we can have Rayleigh, 
(29) Rician, or log-normal channels. Consider first the Rician 
channel, since it is Rayleigh with an added specular compo- 
nent. If we pass transmitted signal (31) through a Rician 
channel we get 


exp (ocg - 0!jX); x>T 


fix) = 


2 ’ 


x<T 


(30) 


r(f) = \/2A sin (co^f + /n(f)] + \/2 n^(f) cos [co^f + m(f)j 
+ \/Tn (t) sin |o) f + m(f)j + n(f) (33) 


and 


a^j = 4.4514 a, = 5.7230. 


where and are zero mean Gaussian random processes 
each with variance . Then the received signal amplitude a(f) 
has the Rician probability density function 


T is the Viterbi decoder node synchronization threshold, „ , , , I a 4\ 

where P(«) = ^ ®xp |-(a" ^A^)l2o^\l^ 


In 2 + a 
T> 

«i 


Define 


and with equality for perfect node synchronization. In this 
paper we have assumed perfect node synchronization. The 
effect of node sync on RS decoding can be found using (25) 
through (30) with the model for n given in (1); results are illus- 
trated in Fig. 9. Further detail is given in Ref. 9. 

For the effect of a noisy carrier reference on RS/Viterbi, 
we also recommend Ref. 5 to the interested reader. 


VII. Effects of Slow Fading on the RS/Viterbi 
Decoder Performance 

Usually the RF telemetry signal can be expressed as 

5(f) = y/2A sin [co^f + /«(f)] (31) 


y 


2 


A ^ 

2or^ 


(35) 


as the ratio of the specular power to the fading power. 
Suppose the received power is P where 


P = A'^ + 2a^ 


(36) 


Let’s normalize a(f) as 


yit) 


g(0 

n/p 


(37) 
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Then 

Piy) = 2;^ (1 + 7^) exp|-(l + - 7 ^} 

/o [2 j'n/7'(1+7^)1 (38) 

with 

7" = 1 (39) 

For the Rayleigh channel we don’t have specular components, 
which means ^ = 0 or 7 ^ =0. Then pdf for y is 

pO) = 2y exp (-y^) y>Q (40) 

Figure 10 shows p(y) given in (38) and (40). For Log-normal 
channels the received signal amplitude is of the form 

a(t) = A (41) 

where x(0 is Gaussian random process with variance a\ Pdf of 
normalized a{t) is 


Using these averages in (2), (5), and (18) we can get 
performance of concatenated coded system in the presence of 
a slow fading channel. The corresponding performance curves 
are shown in Figs. 1 1 and 12. 

Vill. Combined Effect of Phase Jitter and 
Siow Fading 

With a similar argument as before, we can find the 
combined effect of phase jitter and slow fading by finding the 
average of rr and Pj, (Viterbi) over the phase error <t> and 
amplitude fading factor as 

7T = 7r(y,0) = J* j y^ cos^ pWy) d<t>dy 

(46) 

where (for details of effects of a fading channel on PLL see 
Ref. 10). 

p^4,\y) = (47) 

2rrl^(p(y^)) 


Piy) = — exp {-(Iny + I2a\} y > 0 (42) 


pO^) 



N^B^{y^)T{y^) 


(48) 


Therefore in all cases we can assume the received signal is 

r(0 = y/Wy cos + m{t) + 0(t)] + n{t) (43) 

where pdf ofy is given by (38), (40) or (42). 

If the signal amplitude changes very slowly in comparison 
to the bit rate, y will remain constant over a large number of 
bits. On the other hand, 16 levels of interleaving makes y 
affect RS symbols independently. Therefore, we have 

n = iT(y) = J y^j p.(y)dy (44) 

and 

P^iViterbi) = (Viterbily'^) = j y^ p(y)dy 

( 45 ) 


and similarly 

( Viterbi) = J" J y 2 ^^^2 ^ 

(49) 

Using these in our performance formulas, we get the desired 
results. 

IX. Conclusion 

Key parameters for characterizing the performance of 
concatenated Reed-Solomon/Viterbi coding have been con- 
sidered. Simple derivations of close form expressions for a 
number of error probabilities are presented; these include RS 
codeword, RS information block, RS symbol, and RS bit error 
probabilities. A functional model for the RS decoder input 
symbol error probability is found which enables us to carry 
out numerical computations of the above-mentioned error 
probabilities. In addition, the effects of noisy carrier reference 
and slow fading on the RS/Viterbi decoding performance are 
determined. 
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Fig. 1. Concatenated Reed-Solomon/Viterbi coding system diagram 
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CODEWORD SIZE = N = M- I 


Fig. 2. Basic Reed-Solomon code structure 


DMC 



Fig. 3. Equivalent M-ary discrete memoryless channel (DMC) of Fig. 1 
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Fig. 8. Effects of noisy carrier reference on symbol error 
probability 77 


Fig. 7. P„ f(RS), Pg^and P,y performance curves 
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BIT ERROR PROBABILITY (OUTER CODE) P, (RS) 
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Fig. 9. Effects of noisy carrier reference on concatenated Fig. 10. Rician ( 7 ^ >0) and Rayieigh = O)iprobabiiity density 

RS/Viterbi bit error rate (RS ) functions 
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BIT SNR (OUTER CODE), dB BIT SNR (OUTER CODE), dB 

Fig. 11. Pj,(fiS) vs bit SNR for Rician (7*>0) and Rayleigh (y* = 0) Fig. 12. (RS ) vs bit SNR for log-normal fading channels 

fading channels 
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Cn} 


Here we would like to find a simple upperbound on P^. (xq -»■ x„) using the Chernoff 
bound technique. Note for X > 0 we have 


/V{Xj,^x^} = Pr{W + Z<t} = Pr{t- W - Z > 0} < 




(A-1) 


But 


w. 

E{e~^'^} = E(e = 

1 = 1 


- 1 


m- \ \ m- 




(A-2) 


Similarly 


N-n 


E[e-^} = E{e '=> > = [(1 - rr) + Tre"^]^-" 


(A-3) 


Therefore, 


Pr {Xg -> x^ } < min e 




\>o 


— ^ + (l-— ^ 
m - 1 I m - 


[(1 - 7r) + 7re-^]"-" 


(A-4) 
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